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1 Introduction 



For < e < 2, we consider the operator 



W)(*):=4(-»W§Q 




(i) 



initially defined on all C 2 periodic functions / G L 2 {— ir, 7r); the exact domain 
is given by taking the closure of the operator defined on the above functions. 
That such a closure exists will be shown later. In a recent paper |lj Benilov, 
O'Brien and Sazonov showed that the equation 



approximates the evolution of a liquid film inside a rotating horizontal cylin- 
der. 

We shall show that the eigenvalue problem 



has only real eigenvalues, which were shown to exist by Davies in [3]. In the 
same paper, he showed that the spectrum of —iH is equal to the set of its 
eigenvalues, so this implies that the spectrum is real. This was conjectured 
in pQ , and Chugunova and Pelinovsky proved in [2] that all but finitely many 
eigenvalues are real, and gave numerical evidence that all eigenvalues are 
real. Our approach is to show that it is sufficient to consider H acting on 
the Hardy space H 2 (—tt,tt) and analytically continue any solution of (j3J) to 
the unit disc, where the corresponding ODE (jSJ) is now self-adjoint on [0, 1] 



-iHf = Xf 



(3) 




with regular singularities at the end points. In order to do this we make use 
of a bound on the Fourier coefficients proved by Davies in [3J. 
As in [3], by expanding / G L 2 (— tt, tt) in the form 

V ^7T neZ 

one may rewrite the eigenvalue problem in the form Av = Xv, where A = —iH 
is given by 

(Av) n = ^n(n - l)u n _i - |n(n + l)v n+ i + nu n . 

Here we have identified / 2 (Z) and L 2 (— 7r,7r) using the Fourier transform 
T : Z 2 (Z) -> L 2 (-tt, tt). We have 

and define Dom (H) = !F(Dom(A)). 

The (unbounded) tridiagonal matrix A is of the form 



A 



where A_ acts in / 2 (Z_), the central acts in C and A + acts in Z 2 (Z + ). We 
assume that A + has its natural maximal domain 

V = {ve Z 2 (Z+) : A + v e / 2 (Z + )}. 

Davies has shown in [3J that A + is closed and that T> is the closure under 
the graph norm of A + of the subspace consisting of those v G ^ 2 (Z + ) that 
have finite support. Let r be the natural identification between / 2 (Z + ) and 
Z 2 (Z_); then Dom (A) = r(T>) © C © V and r induces a unitary equivalence 
between A + and A-. Therefore, in order to prove that all eigenvalues of 
A are real, we only need to prove that all eigenvalues of A + are real. The 
Fourier transform identifies / 2 (Z + ) with {/ G H 2 (-n,ir) : Jl 7T f(e)de = 0}. 

Let H be the restriction of H to C 2 er ([— n, tt]), which is clearly a subspace 
of Dom (H). We now show that H is the closure of Hq. 

Proposition 1.1 Where H and H are as above, H is the closure of Hq. 

Proof It follows from Davies' result on the domain of A + that the trigono- 
metric polynomials are dense in Dom (H) with respect to the graph norm. 
Since the trigonometric polynomials are contained in C 2 er ([— tt, tt]), this space 
is also dense in Dom(H), which is closed in graph norm since Dom (A) is. 

I 
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2 Reality Of The Eigenvalues 

If A + v = Xv, then v is a solution of the recurrence relation 

Un + l)(n + 2)t» n+2 + (A - n - l)v n+1 - S -n{n + l)v n = (4) 

satisfying the initial condition ev 2 = (1 — X)v\. We shall study the generating 
function, Y?£=i v k zk °f { v k}- 

Lemma 2.1 Let v G l 2 (Z + ) be such that A + v = Xv. Then the function 
u(z) := SfeLi VkZ k , defined for \z\ < 1, satisfies the differential equation 

„ z + l/e , 2\/e , , 

u" - 2- ±— ^ r-u' + — -u = 0. 5 

(l-z)(l + z) z(l-z)(l + z) v 7 

Proof The constant term in 

^(l - z)(l + z)u" - 2(z + l/e)zu' + (2A/e)« 
is clearly 0, and the coefficient of z is 

2v 2 - 2vi/e + 2A^i/e = 2 ^ 2 - - ~ A ^ = 0. 

The coefficient of z n is 

2 2A 
n(n + l)u n+ i - (n - l)(n - 2)w n _i nv n - 2(n - l)u„_i H v n 

£ £ 

= - |n(n + l)v n +i + (A - n)w n - |n(n - l)v n -i = 
for n > 2. I 

From here on we assume that A G C is an eigenvalue of the operator A + 
and v is a corresponding non-zero eigenvector. In Davies proved that 
there exist constants 6, m such that 

IMU<W"iMI 2 , (6) 

where c = 1 + 1/e and | j"^ I loo c := sup{|w n |n c : 1 < n < 00}. Since e > 0, this 
implies that v G / 1 (Z + ), and hence that J2^i v kZ k is absolutely convergent 
for \z\ < 1. The equation (jSJ) can be written in the form of Heun's equation 

U + - + + U + — — rU = (7) 

\z z—1 z — aj z\z — l)(z — a) 
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with a = 1, (3 = 0, 7 = 0, 5 = 1 + 1/e, e = 1 — 1/e, a = — 1 and = 2A/e. 
This is a Fuchsian equation with four regular singular points, at 0, 1, — 1 
and oo, with {0, 1 — 7}, {0, 1 — 5}, {0, 1 — e} and as the roots of 

the corresponding indicial equations (for the Frobenius series at each regular 
singular point). For background information on Heun's equation, see [I]. 

Lemma 2.2 Suppose that < e < 2, 1/f ^ Z, A 6 C is an eigenvalue of A + 
and v G ^ 2 (Z + ) is a corresponding non-zero eigenvector. Then there exists a 
solution u of (0) which is analytic in an open set containing [0, 1] and such 
that u(z) = J2kL\ v kZ k for all z such that \z\ < 1. 

Proof Define u(z) = YX=i v k zk on {z G C : \z\ < 1}. Let u\ be the solution 
of (jTJ) with exponent about 1 and be the solution with exponent —1/e 
about 1. Let U be the intersection of the open discs of unit radius about 
and 1. The space of solutions of ((7|) in U is two-dimensional, and u, U\ and 
U2 lie in this space. Hence there exist constants a, b such that u = au\ + bu<i 
in U. Since v G / 1 (Z + ), u(z) converges to a finite limit as z — > 1 in U. Also 
U\{1) is finite, but U2{z) — > 00 as z — > 1 in U. Therefore we must have b = 
and u = aui in U. Let W be the union of the open discs of unit radius about 
and 1. We now extend u to all of W by u(z) = aui(z) on the open disc 
of unit radius about 1. Now u is an analytic solution of (171) on W, which is 
an open set containing [0, 1], such that u(z) = J2T=i v k zk f° r ai l z such that 
\z\ < 1. I 

Theorem 2.3 Suppose that < e < 2, 1/e ^ Z and A G C is an eigenvalue 
of A+. Then A G R. 

Proof Let u G ^ 2 (Z + ) be a non-zero eigenvector corresponding to A. Let u 
be as in Lemma 12721 and put fi = 2X/e. Following [5], the equation (JTJ can 
be written as 

— (pu 1 )' + qu = [muu (8) 
on the complex plane cut along [1, 00) and (—00, 0], where 

p(z) = z" 1 {1 - z) s (z - af 
q\z) = a(3z^(l- z) 5 - l (z - a) e ~ 
w( z ) = z T-i(i _ z y-i( z - a) 6 " 1 

We now restrict u to [0, 1]. It is clear that u G C°°([0, 1]) and 

— {pu')' = fiWU 



= (l-z) 1+1 / £ (z + l) 1 - 1 / £ 
= 

= z-Hl-zf/Hz + l)- 1 /*. 
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on (0, 1). Note that p > on [0, 1) and w > on (0, 1) with p(x), w(x) — > 
as x — > 1 from below. Since u has a zero of order 1 at and w has a pole of 
order 1 at 0, wu G C([0, 1]). Therefore |-u| 2 u> G L 1 (0, 1). Now 



1 f rl-l/n 



it / lit (a;) I w(x)dx = lim < — / (pu ) (x)u(x)dx 

JO n -*°° I Jl/n 



lim 



lim 

n— >oo 



1-1/n rl-l/n 



p(x)u' (x)u(x) + / p(xW(xW(a:)d£ > 

Jl / n Jl/n J 
, i -i 1— 1/n rl— 1/n 

^ — -u'fxWa;) f — / (jw'Vfa; Wandas 

L J J 1/n Ji/ n 







(pu')'(x)w(x)d£ 



= /i / |w(a;)| 2 w(a;)da;. 

JO 

Since w is a non-zero solution of ([H]) and w > a.e. we have /i = /I and hence 
(iGR. Since A = we also have A G R. | 
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